Simple model of anisotropic pairing with repulsive interactions. 
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A simple tight binding model with repulsive interactions is studied. The inclusion of more than 
one orbital per site leads to assisted hopping effects, and, when the orbitals involved have different 
symmetries, to an anisotropic superconducting phase. Superconductivity exists for all fillings, and 
for all values of the on site repulsion. 
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Introduction. The existence of superconductivity in 
models of strongly correlated electrons with only re- 
pulsive interactions has attracted a great deal of at- 
tention in recent times. It was established, after the 
formulation of the BCS theory, that the metallic state 
is unstable towards anisotropic superconductivity, due 
to the angular dependence of the dielectric constantEl. 
This Kohn-Luttinger instability is areatly enhanced when 
the Fermi surface in anisotropicBBa. By using RPA 
or model dielectric functions, it can also be shown 
that isotropic Feppi surfaces give rise to anisotropic 
superconductivityEl. Varied numerical evidence suggests 
that models with purely repulsiMe,|ijBteractions can lead 
to anisotropic superconductivitj^Blatl. 

An alternative scheme which leads to a superconduct- 
ing ground state starting frnm .models with repulsive in- 
teractions was proposed intHlEil. This model includes an 
assisted hopping term which arises naturally when con- 
sidering many orbitals at each lattice site. This assisted 
hopping term strongbr.fay.Qi-s the existence of a supercon- 
ducting ground stateEJOtJ. In its standard version, this 
model leads to isotropic superconductivity for small hole 
concentrations, and for a moderate value of the onsite 



repulsion. In this work, we present an extension of this 
model which has an anisotropic superconducting ground 
state for all dopings and arbitrary values of the onsite 
repulsion. 

The model. We study the simplest model with the fea- 
tures discussed above. We postpone to the conclusions 
the analysis of possible generalizations. Following the 
scheme inliJ, we assume an atom with two orbitals per 
site on a slightly distorted square lattice. The distor- 
tion makes the lattice orthorhombic. We assume that 
this distortion is small, and we will treat its effects as 
a perturbation. Using the tetragonal symmetry of the 
undistorted lattice, we take the lowest lying orbital to be 
d^2_y2 , and the second orbital to be s or d^^2_j.2 . We de- 
fine A = Cs — ea;2_y2 as the difference between the energy 
of these levels. Without loss of geenrality, we set — 0. 
We assume that there is only repulsion between electrons 
in the dx2_y2 orbital, U. Hopping can take place only be- 
tween nearest neighbors, with amplitudes tss, tad and tdd- 
Finally, crystal field effects associated to the orthorhom- 
bic distortion induce an hybridization between the s and 
dx2_y2 in the same site, Vcf- The hamiltonian is: 
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For simplicity, we neglect here possible differences be- 
tween the hoppinggLin the two directions due to the asym- 
metry of the latticeE3, which is included through the crys- 
tal field potential Vcf only. Note that, due to the differ- 
ent symmetries of the two orbitals, the hopping between 
an s and a d orbital has opposite sign along the two axes 
of the lattice (see fig.||l|). 

We now assume that Vcf ^ U, \tss\, \tsd\, \tdd\ ^ A. 
The lowest lying eigenstates of Tiion for each occupancy 
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with energies: 
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Using eq.(g) we define an effective hopping which de- 
pends on tlic occupancy of the site. We find: 



tunn 

{^^i\{^lj\'Htunn\2i)\0j) ~ tdd^tsd 
(lcri|(2j|Htunn||2i)lo-j) ~ tdd^'^tsd 



VCF 

A 

fVcF 



CF 



\ A A-U 

VcF 

A-U 



(4) 



FIG. 1: Sketch of the hopping terms between the two orbitals 
in the unit cell. 



We can now write an effective hamiltonian: 



i aij 
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where we have shifted the origin of energies by —V^p/A, 
see eq.(||), and: 

~ ^ mY£L 
St « 2Ud^j- 

U^U-2^ (6) 

Note that, as we are expanding to first order in Vcf/A, 
we can neglect normalization terms in the definition of 
the electron operators c^i. For the same reason, 6t has 
the same absolute value along the two axes of the lattice. 
The symmetries of the orbitals involved imply that, to 
first order, there are not next nearest neighbor assisted 



hopping terms. 

Superconducting solution. The energy bands of the ef- 
fective hamiltonian, eq.(^, are: 

ek^ky = {i+ + ^Stj cos(fcj;) + (t^ - ^Sij cos(fcy) (7) 

Where n is the number of electrons in the unit cell. 
The pairing interaction isll3: 

Vk^kyk'^k'y = U+Si [cos(fcj;) - cos{ky) + cos(fc^) - cos(fc,y)] 

(8) 

The superconducting gap must be of the form: 

Asc k^ky = a + b [cos{k^) ~ cos(A;j,)] (9) 

Following ref.lli, at the transition temperature these co- 
efficients satisfy: 



a = a Ul0-6i{l^- ly) +b U{I^ - ly) + 5t{L^^ - 2I^y + lyy) 

b = aSilo - b6i{h - ly) 
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where: 



cos{k^)n{ek^kjT) 

k^ky Cfexfcy — 



and similar expressions for ly , Ixy and /yy . The function 
n{uj/T) is the Fermi-Dirac distribution, and e^? is the 
Fermi energy. 



(11) From eqs.(lO), the critical temperature is given by: 



= 1 + f//0 - 25i{I, - ly) + SP [(/, - ly)^ ~ h{h^ - 2I,y + lyy)] 



(12) 



The integrals in eqs.(|ll|) diverge as — log[{N{eF)T] at 
low temperatures, where N{eF) ~ is the density of 
states at the Fermi level, and i = (t+ = i-)/2. As we 
are assumming a weak orthorhombic distortion, ^ ly 



Hence, eq.(|12|) always has a solution at low temperatures. 
The critical temperature is given, approximately, by: 



Tc ~ Cite 



(13) 

where we assume tdd ~ tsd ~ t. c\ and ci are numer- 
ical constants. The expression in eq.(|l3|) implies that 
the superconductivity is supressed as A ^ cx). It is en- 
hanced by the the lattice asymmetry, described by VcFi 
and by the existence of a finite ?7, as the assisted hop- 
ping term requires the presence of electron-electron in- 
teractions. Near half filling, the van Hove singularity in 
the density aLjStates implies cx te~°^, where A ~ 
A'*/(iV(^^C/)BBE3. The expression in eq.(|l|) ceases to 
be valid near the band edges, as I^x — 21, 
can be shown that, in this limit: 
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where n is the number of carriers per unit cell. Then: 



We can also calculate: 
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so that a/5 —> as A becomes the largest energy in the 
problem. In this limit, the gap, eq.(^), will have dx2-y2 
symmetry. 

It is interesting to compare the value in eq.dl^) 
with the critical temperature deduced from the Kohn- 
Luttinger analysis for the Hubbard model without as- 
sisted hopping terms, T^kl- The effective coupling con- 
stant arises from the screened interaction, and, to low- 
est order in U, it goes as U'^N{eF)- Hence, Tckl ^ 
dite"'*2/(^^('^)l'. Thus, within the perturbative ap- 
proach used here, we find Tc ^ Tckl, although Tc also 
depends on the strength of the orthorhombic distortion. 



Phase diagram. So far, we have only considered the 
superconducting instability. Near half filling, it is well 
known that the hamiltonian in eq.(|l|) has nesting prop- 
erties, and antiferromagnetism is favored. The corre- 
spondiM Neel temperature is proportional to the gap, 
so thatO: 



T-i 



N 



(17) 



Assuming that Vcf ^ A, we find that, near half fiUing, 
Tc <C Tat, and the system will be antiferromagnetic. This 
phase dissappears at densities, measured from half filling, 
n ^ Q-i'^i/ \/VN(tF) _ sketch of the expected phase dia- 
gram is shown in Fig.|^. Th method used here cannot 
be used to elucidate the nature of the transition from 
the antiferromagnetic phase at half filling to the super- 
conducting phase. Hartree-Fock studies of the Hubbard 
model suggest that an orthorhombic distortion favors the 
formation of_stripeai3. Alternatively, phase separation is 
also possiblell3. Higher order terms in Vcf/ U/A, not 
considered herei—will break the electron-hole symmetry 

m. It is 



shown in Fig. 



is interesting to note that the sin- 

cy 



gie band Hubbard model can show an intrinsic teui 
towards an orthorhombic distortion near half fiUiu] 

Conclusions. We have studied a variation of the as- 
sisted hopping model, which includes tiui juffects of other 
atomic levels in Hubbard hamiltoniansEJ'Eil (see alsocj). 
We have avoided the difficulties associated to intermedi- 
ate coupling situations by restricting the study to a well 
defined weak coupling case, where standard perturbation 
theory is reliable (for |a,similar approach to other assisted 
hopping problem, seeu). 

The induced assisted hopping terms lead to d-wave su- 
perconductivity. The arguments used here show that the 
symmetry of the order parameter, for models with dif- 
ferent atomic orbitals, is associated to the relative sym- 
metry of these orbitals. When the most relevant orbitals 
are of the same symmetry, onlyjS^mave superconductiv- 
ity is possible, in agreement withE3llj. It is interesting to 
note that the symmetry of the order parameter does not 
arise from the anisotropics of the Fermi surface, as for the 
Kohn-Luttinger mechanism. We have shown, that in the 
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FIG. 2: Sketch of expected phase diagram for the model 
described by eq.(Q), see text for details. The present anal- 
ysis is insufficient to characterize the intermediate regions 
between the antiferromagnetic and superconducting phases, 
shown shaded in the figure. 



weak coupling regime, the value of the critical tempera- 
ture arising from the assisted hopping term tends to be 
larger than that due to the Kohn-Luttinger mechanism. 
It is plausible that, in the same manner, the presence of 
assisted hopping terms will enhance the tendency of the 
standard Hubbard model towards superconductivity in 
the intermediate coupling regime. 

Finally, the present results can be extended to other 
systems where assisted hopping terms are likely to arise. 
In quaiitwm dots, where the sign of these terms can be 
randoniSil, they can enhance the tendency towards local 
pairing. 
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